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ABSTRACT 

Context. Characterizing the spatial and velocity structure of molecular clouds is a first step towards a better understanding of interstellar 
I turbulence and its link to star formation. 
f* «) ■ Aims. We present observations and structure analysis results for a large-scale (~ 7.10 deg 2 ) 13 CO J = 2-1 and 12 CO J = 3-2 survey towards the 
\^ | nearby Perseus molecular cloud observed with the KOSMA 3m telescope. 

C") Methods. We study the spatial structure of line-integrated and velocity channel maps, measuring the A- variance as a function of size scale. We 
determine the spectral index ft of the corresponding power spectrum and study its variation across the cloud and across the lines. 
Results. We find that the spectra of all CO line-integrated maps of the whole complex show the same index, « 3.1, for scales between about 
0.2 and 3pc, independent of isotopomer and rotational transition. A complementary 2MASS map of optical extinction shows a noticeably 
O smaller index of 2.6. In contrast to the overall region, the CO maps of individual subregions show a significant variation of p. The 12 CO 
3-2 data provide e.g. a spread of indices between 2.9 in L 1455 and 3.5 in NGC 1333. In general, active star forming regions show a larger 



6 

ffi ' power-law exponent. We find that the A-variance spectra of individual velocity channel maps are very sensitive to optical depth effects clearly 
. . indicating self-absorption in the densest regions. When studying the dependence of the channel-map spectra as a function of the velocity 
channel width, the expected systematic increase of the spectral index with channel width is only detected in the blue line wings. This could be 
explained by a filamentary, pillar-like structure which is left at low velocities while the overall molecular gas is swept up by a supernova shock 
wave. 
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1. Introduction 

The structure of the interstellar medium (ISM) is random 
to a large degree, with a complex spatial density distribu- 
tion and velocity field. This is evident in large-scale sur- 
veys of spectral line tracerS j such as CO, m a de for parts 
of th e Galactic plane (e.g. iHever et all Il998t Simon et ail 
20011) and individual molecula r clouds (e.g. lBallv et alJll987t 
Ungerechts & Thaddeus 1981 . The spatial structure of the 
emission has been quantified in terms of its powe r spectrum 
dScalolll987t fstutzki et aljfl998t lElmegreer]ll999h . When fit- 
ting the azimuthally averaged power spectrum with a power 
law, the slope of the power law j3 provides information on 
the relative amount of structure at the linear scales resolved 
in the image. A pure power law is expected for structure 
on the linear scales of a self-similar turbulent energy cas- 
cade. Deviations from a power law are expected at scales 
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where physical processes insert energy in the turbulence cas- 
cade (outflows, supernovae, super-bubbles , galactic rota tion) 
and a t scales of turbulence dissipation ( Elm egreen & Scale! 
12004 . Thus, a study of the scaling behavior of the cloud 
structure and the velocity field as traced by the power spec- 
trum of observed spectral line maps can help to constrain the 
turbulent energy cascade in the ISM. A number of power- 
spectrum studies have b een carried out for the atomic med ium 
using Hi observations (Crovisier & Dickev 1983; Greenll993l 
I Stanimirovic et all 1 19991 iDickev et alJl200ll lElmegreen et alJ 
l200lh . They find spectral indices between abou t 2.7 and 3.7 for 
the line-integrated maps ( Fal garone et al 120041) . 

The A-variance is an alternative method to determine the 
index of the power spectrum of isotropic images JStutzki et alJ 
1998). In contrast to the power spectrum, the A-variance can 
be computed in the spatial domain. It allows for a better sep- 
aration of the intrinsic cloud structure from contributions re- 
sulting from the finite signal-to-noise in the data and the tele- 
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scope beam. In addition, problems related to the d i screte sam- 
pling of the data can be avoided. iBensch et al .] d200 111 pre- 
sented a A-variance analysis for CO maps of the Polaris Flare 
and six other molecular clouds. The index f3 is close to 3 for 
most clouds, but it steepens in the Polaris Flare from 2.5 to 
3.3 for maps with a linear resolution increasing from k. 1 pc 
to ^ 0.1 pc. lOssenkoo f et al.l (|2£)01) used the A-variance to 
study numerical models of self-gravitating, supersonically tur- 
bulent clouds and compared the results to observations. They 
showed that the A-variance traces deviations from an inertial 
scaling behav ior at the scales of driving and dissipation (see 
also Ossenkopf & Mac Lowl2002l) . 

Constraints on the velocity structure of the ISM can be ob- 
tained from an analysis of individual channel maps of a spectral 
line cube, and by comparin g the results with those ob tained 
for the line-integrated maps. Stanimi rovic & Lazari an (2001) 
have analyzed H i observations of the Small Magellanic Cloud 
in this way, noting a systematic variation of the measured index 
j6 with the width of th e velocity-channels, thus confi rming theo- 
retical predictions by Lazarian & Pogosvan (2000). Using the 
Boston University/Five College Radio Astronomy Obse r vatory 
(BU/FCRAO) Galactic ring survey, IBensch & Simonl feOOlh 
found that the index of the channel maps is smaller than that 
of the line-integrated maps. It showed a significant variation as 
a function of the channel velocity. But [} did not increase with 
in creasing velocity channel widt h, in contrast to the predictions 
byE azarian & Pogosvan (2000). Apart from this first attempt, 
no systematic structure analysis of velocity-width dependent 
channel maps of molecular lines has been performed yet. 

We present a new CO survey of the Perseus molecular 
cloud complex and analyze the cloud structure in the line in- 
tegrated and the channel maps. The observations of the 12 CO 
3-2 and 13 CO 2-1 transitions were made with the 3m telescope 
of the Kolner Observatorium fur Sub-Millimeter Astronomie 
(KOSMA) and cover the entire Perseus molecular cloud com- 
plex (7.1 square degrees). The Perseus cloud is one of the best 
examples of a nearby star-forming regio n. The cloud is re- 
lated to the Perseus OB 2 association JBachiller & Cernicharol 
1 19861 llJngerechts & Thaddeus 1987), which is located in the 
area mostly free of CO emis sions northeast of the cloud 
JlJngerechts & Thaddeusl fl987t) . It includes a region where 
intermediate-mass stars form (NGC 1333), a young open clus- 
ter (IC 348), and a dozen dense cloud cores with low levels 
of star-formation activity (L 1448, L 1455, B 1, B 1 EAST, B 3 
and B 5). 91 protostars and pre-stellar cores have been identi- 
fied in a 3 square degree survey of the dust continuum at 850 
and 450 i m\ made with the James Clerk Maxwell Telescope, 
JCMT llHatchell et all2005h . 

Imaging observations of the Perseus complex in molecular 
cloud tracers exhibit a wealth of substructure, such as cores, 
shells, filaments, outflows, jets, and a large-sca le velocity gra- 
dient JPadoan etalJ Fl999) . Pado an et alJ i 19991) compared the 
structure traced by 13 CO 1-0 observations to synthetic spectra 
and find that the motions in the cloud must be super- Alfvenic , 
with the exception of the B l core, where Good man et al.l fl989) 
and Crutcher et al Ji \l 9931) detected a strong magnetic field. 
IPadoan et alJ(l2003al) find that the structure function of the line- 
integrated 13 CO 1-0 map follows a power law for linear scales 



between 0.3-3 pc, and Pad oan et alJ ll2003bl) compared the ve- 
locity structure of Perseus to MHD simulations. 

A complete census of the stellar content of nearby 
(^ 350 pc) molecular clouds (Perseus, Serpens, Ophiuchus, 
Camaeleon, and Lupus) is currentl y obtained by the Spitzer 
legacy project "Cores to Disks" (c2d. lEvans et alJ2003l) . Large- 
scale 12 CO, 13 CO 1-0 and A v maps of the nor thern clouds wer e 
recently obtained by the COMPLETE team JGoodmanl20 04). 
The KOSMA survey of Perseus in higher CO transitions traces 
the warmer and denser gas due to the elevated critical densities 
and excitation energies (~ 10 5 cirT 3 and 33.2 K for CO 3-2) 
relative to the J = 1-0 transition. Moreover, 12 CO is largely 
optically thic k, while 13 CO, being a factor ~ 65 less abundant 
dLanger&Penziaslll990l) . is often optically thin, thus tracing 
column densities. We plan to extend this work to the other 
nearby clouds using the KOSMA and NANTEN2 1 observa- 
tories. 

The KOSMA observations are described in Sect. 2. The 
general properties of the CO data sets are discussed in Sect. 
3. Sect. 4 presents the results of the A-variance analysis. The 
discussion of the results and a summary are given in Section 5 
and 6, respectively. 

2. Observations 

We mapped the Perseus region simultaneously in 13 CO 2-1 
and 12 CO 3-2 using the KOSMA 3m submillimeter telescope 
on Gornergrat , Switzerlan d, equipped with a dual-channel SIS 
receiver ( Graf et al. | l l998l) and acousto optical spectrometers 
( Sch ieder et al.lfl989l) . Main beam efficiencies and half power 
beamwidths (HPBWs) are 68%, 130" at 220 GHz and 70%, 
82" at 345 GHz. The HPBWs correspond to linear resolu- 
tions of 0. 22 pc and 0.14pc, where we adopted a distance 
of 350pc JBorgman& Blaauwlll964l: iHerbig & Joneslll983t 
iBachiller & Ce rnicharo 1986fKAll temperatures quoted in this 
paper are given on the main beam temperature scale. 

For the observations, we divided the ~7.10 deg 2 re- 
gion in 10' x 10' fields. Each field was mapped using the 
position-switched o n-the-fly (OTF) mode JKramer et al.l ll999: 
Beut her et alJ l2000) with a sampling of 30". Three emission- 
free off positions were selected from the I3 CO 1-0 FCRAO 
map. The pointing was accurate to within 10". The accuracy 
of the absolute intensity calibration is better than 15%, de- 
termined with frequent observations of reference sources. The 
channel spacing Av c h and the corresponding average baseline 
noise rms of the spectra is 0.22km s -1 , 0.48 K for 13 CO 2-1 
and 0.29 km s~', 1.02 K for 12 CO 3-2. The observations were 
taken from February to December 2004. 

3. Data Sets 

3.1. Integrated intensity maps 

The maps of velocity integrated 13 CO 2-1 and I2 CO 3-2 emis- 
sion (Figs. [33 show the Perseus region, viz. the well known 
string of molecular clouds running over ~ 30 pc projected dis- 
tance from NGC 1333 and L 1455 in the west to B 1, B 1 East, 

1 http://www.ph 1 .uni-koeln.de/nanten2 
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and B 3 in the center, and to IC 348 and B ,5 in the east (cf . 
iBachiUer & Cernicharolll986t lUneerechts & Thaddeudll987t> . 
Generally, there is a good correlation between 12 CO 3-2 and 
13 CO 2-1 integrated intensities. 

In the next section, we compare the statistical properties of 
the structure seen in the entire Perseus map with the structure 
seen in individual regions. For this, we defined seven boxes of 
50' x 50' which roughly coincide with the known molecular 
clouds (cf. Fig.^. 

Figure [2] shows an overlay of integrat ed 13 CO 2-1 in- 
tensi ties and a m ap of optical extinctions 
lAlves et al]Eo 05). at 2.5' and 5' resolution, respectively. The 
13 CO map covers all regions above 7 mag and ~ 70% of the 
regions above 3 mag. A linear least squares fit to a plot of A,, 
vs. 13 CO 2-1 results in a correlation coefficient of 0.76. The 
region mapped in 13 CO has a mass of 1.7 x 10 4 M Q using the 
Ay data and the can onical conversion f actor [H2]/[Ay] = 9.36 
1 llBohlin et alii 9781). 



x 10 20 cm 2 mag 



3.2. Velocity structure 

Maps of 13 CO 2-1 emission integrated over small velocity in- 
tervals (Fig.[3} illustrate the filamentary structure of the Perseus 
clouds. The channel maps show the well-known velocity gradi- 
ent between the western sources, e.g. NGC 1333 at ~ 7 kms -1 , 
and the eastern sources, e.g. IC 348 at ~ 9 km s _I . The channel 
map integrated between 5 and 6kms~ 1 exhibits two filaments 
originating at L 1455, one runs north to NGC 1333, the second 
runs north-east to B 1 . We will discuss the structural properties 
of individual velocity channel maps in the next sections. 

To study the statistics of the velocity field we start with 
the distribution of the line widths across the map. Since many 
spectra show deviations from a Gaussian line shape, we use the 
equivalent line width Av eq = J Tdv/T ped k as a measure of the 
line dispersion along individual lines of sight. Figure |4] shows 
the mean equivalent line widths and their scatter for the seven 
sub-regions shown in Figure^ 

The mean I2 CO widths vary significantly between 
2.2kms~ 1 in the quiescent dark cloud L 1455 and 3.8kms _I 
in the active star forming region NGC 1333, while the rms is 
~ 0.7 km s _1 . In contrast, the I3 CO widths are smaller and show 
only a weak trend around ~ 2 kms 1 . 

Several positions in L 1455, but also in e.g. IC 348, show 
small line widths of ~1 kms -1 , only a factor of — 8 — 1 1 larger 
than the CO thermal line width, which is « 0.16 kms 1 for a 
kinetic temperature of 10 K as was fou nd for the bulk of the gas 
in Perseus by Bachil ler"& Cernicharol {l^ 
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Fig. 4. Mean and rms of the equivalent line widths Av eq of the 
12 CO 3-2 and 13 CO 2-1 spectra for the observed positions of 
the seven 50' x 50' sub regions (Fig.^. The dashed line delin- 
eates equal widths in 12 CO and 13 CO. 



All Sky Survey) by the COMPLETE team JCioodmanl 12004: 
lAlveset.all200.il 

4.1. The A-variance analysis 

The A-variance analysis was introduced bv lStutzki etaDdl998l) 
as a means to quantify the relative amount of structural varia- 
tion at a particular scale in a two-dimensional map or a three- 
dimensional data set. The A-variance is defined as the variance 
of an image s(r) convolved with a normalized spherically sym- 
metric wavelet O of size L 



o-\ = ([s(r) * © L (r)] 2 > r , 



(1) 



where the asterix denotes the spatial convolution dStutzki et alJ 
1998). For structures characterized by a power-law spectrum, 
P(\k\) oc \k\~P, the A-variance follows as well a power law, with 
the exponent d& — f3 - 2 in the range < ft < 6. 

Unfortunately, the A-variance spectrum of any observed 
data set does not only reflect the spectral index ft of the astro- 
physical structure, but it is al so affected by rad iometric noise 
and the finite telescope beam jBensch et all20oil) . Both effects 
change the spectrum at small scales. When we ignore the small 
contribution from beam blurring, we can write the A-variance 



4. Structure analysis 

In this section, we statistically quantify the spatial structure ob- 
served in the maps, both for the overall structure and for the 
structure of individual regions within the Perseus molecular 
cloud. We measure the spectral index of the power spectrum 
using the A-variance analysis, a wavelet convolution technique. 
We analyze our new CO data and compare the results with an 
equivalent analysis of the FCRAO 12 CO 1-0, I3 CO 1-0 maps 
and the Ay Perseus map obtained from 2MASS (Two Micron 



cr 2 (L) * aiL< 3 - 2 + a 2 L 



(2) 



(Eq. (10) from lBensch etai]l200ll) . For white noise, j8 no j se = 
0, so that c/ no i se = -2. In the KOSMA data we noticed that 
the noise does not follow a pure white noise behaviour, but it 
is "colored" due to artifacts from instrumental drifts, baseline 
ripples, OTF stripes etc. This has to be taken into account when 
deriving the cloud spectral index ft from the A-variance spectra. 

Thus we measured the spectral index of the colored noise 
<^noise by analyzing maps created from velocity channels which 



4 



Sun et al.: A KOSMA 7 deg 2 CO 2-1 & 3-2 survey of the Perseus cloud 



arcmin 

130 IBS 9(1 



1 1 1 


i | i | i | i | i | i | i 


i 1 i 1 i 1 i 1 i 


i | i | i | i 


1 1 1 1 1 


1 i _ 






10po 














E fcm/B 




~ In * 






B ^ , i , 


i i i 1 i_ 






HHnlfl / B 3 




1 


ID 1 




B 5 










(fff^ ' / 




NGC 1333 




— 




B 










BH-_ " ' ^HH 




































pa 






IBSsfl 




































B 1 EAST 


















- ResoL 150" 




















144B ~~ 


— • 






L 1455 HBgl 






; 1 , 1 , 1 , 


1,1,1,1,1,1,1 


,1,1,1,1, 


1 , 1 , 1 i, 1 







1ST 



iar co 1 jef ae" 2j" 

R,A. (1950) 



lflT 



Fig. 1. The Perseus molecular cloud complex. KOSMA maps of integrated intensities of 13 CO 2-1 (colors) and 12 CO 3— 
2 (contours) at 150" resolution. The integration interval is 0-16 km s~'. Colors run from lKkms 1 (~ 1 cr) to 32Kkms _1 . 
Contours range from 6.6Kkms~ 1 (~ 3 cr) to 83 Kkms 1 in steps of 9Kkms _1 . The (0,0) position corresponds to RA=03:26:00, 
DEC=+3 1:10:00 (B1950). Seven sub-regions are marked by dotted square boxes of 50' x 50'. 



do not see any line emission but which cover the same velocity 
width as the actual molecular line maps. The result is shown in 
Fig. |5] We find a nearly constant index c/ no i se ~ —1.5 for all off- 
line channels at scales between about 1 and 6'. At larger lags, 
the noise deviates from the /3 = 0.5 behaviour, but this does not 
affect the structure analysis as the absolute noise contribution 
is negligible there. 

For the FCRAO data and the COMPLETE A v map we have 
no emission-free channels available so that we cannot perform 
an equivalent noise fit there. The A-variance at small lags shows 
however no indications for a deviation from the pure white 
noise behaviour, so that we stick to d no i se = -2 for the fit of 
these data. 

4.2. Integrated intensity maps 

Figure [6] compares the A-variance spectra of the different in- 
tegrated intensity maps for the entire region mapped with 
KOSMA (see Fig. [Q 2 . When corrected for the observational 
noise, the A-variance spectra of all maps follow power laws 
between the linear resolution of the surveys and about 3pc 

2 Note that the area covered by the FCRAO is slightly smaller than 
that observed with KOSMA. 



(Table^. The good agreement of the spectral indices obtained 
from the different CO data is remarkable. They cover only the 
narrow range between 3.03+0.14 and 3.15+0.04. In contrast, 
the extinction data result in a significantly lower index. This 
indicates a more filamentary structure in Ay. When we actually 
compare the Ay map with 13 CO 2-1 data smoothed to the same 
resolution, it is also noticeable by eye that the Ay map looks 
more clumpy or filamentary than the 13 CO map. This indicates 
that 13 CO does not trace all details of the cloud structure, but 
rather measures the more extended, and thus more smoothly 
distributed gas. 

All A-variance spectra show a turnover at about 3 pc. To 
test whether this peak measures the real width of the Perseus 
cloud or whether it is produced by the elongated shape of the 
CO maps, we have repeated the A-variance analysis for the Ay 
data of the entire region shown in Figure |3] In this case we 
find almost the same spectrum below 3 pc, but instead of a 
turnover only a slight decrease of the slope at larger lags. Thus 
we have to conclude that the A-variance spectra of the CO maps 
at scales beyond 3 pc are dominated by edge effects, due to the 
shape of the maps, so that these scales should be excluded from 
the analysis. In Figure|6]we compare only spectra for the same 
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Fig. 2. Overlay of I3 CO 2-1 integrated intensities (contours) with a map of optical extinctions in colors dGoodmanll200i 
lAlves et all20 05). Contours range from 2.7 Kkms 1 (3 <r) to 32 Kkms 1 by 3 Kkms 1 . Colors range from A v = 1 mag to 1 1 mag. 
Resolutions are 2.5' for 13 CO and 5' for Av- A polygon marks the boundary of the I3 CO map. 



region, i.e. the A-variance spectra of the Ay data of the region 
also mapped with KOSMA. 

As it is not guaranteed that the structure of the overall re- 
gion is representative for individual components, we have also 
applied the A-variance analysis to the KOSMA data of the in- 
dividual clouds contained in the seven 50' x 50' subregions 
shown in Fig. ^ The results of the power-law fits to the A- 
variance spectra are listed in Table |2 They differ significantly 
between the individual regions. The active star-forming region 
NGC 1333 shows the highest spectral indices in both transi- 
tions. The low end of the spectral index range is formed by 
the dark cloud L 1455 together with the environment of the 
young cluster IC348. The A-variance spectra of the two ex- 
treme examples NGC 1333 and L 1455 are shown in Fig.[6j>. 
Starting from the same noise values at small scales the spectra 
of the two regions show an increasing difference in the rela- 
tive amount of structure at large scales reflected by the strongly 
deviating spectral indices. Altogether, we find high indices as 
characteristics of large condensations for the regions with ac- 
tive star formation and lower indices quantifying more filamen- 
tary structure for dark clouds, but IC 348 as an exception to this 
rule, showing also a very filamentary structure. 



Table 1. Results of the A-variance analysis of the integrated CO 
maps and the A v data for the region mapped with KOSMA (Fig.0. 



Transition 


Telescope 


resol. 
['] 


Fit Range 
['] 


P 




A v 


2MASS 


5 


5.0-28 


2.55 ± 


0.02 


13 CO 1-0 


FCRAO 


0.77 


0.8-28 


3.09 ± 


0.09 


12 CO 1-0 


FCRAO 


0.77 


0.8-28 


3.08 + 


0.04 


13 CO 2-1 


KOSMA 


2.17 


2.2-28 


3.03 ± 


0.14 


12 CO 3-2 


KOSMA 


1.37 


1.4-28 


3.15 ± 


0.04 



Table 2. Results of the A-variance analysis of the KOSMA data for 
seven 50' x 50' sub-regions of the cloud (Figure^- The spectral indices 
fi were fitted in the size range 2.2-14' for the I3 CO 2-1 and in the size 
range 1.4-14' for the 12 CO 3-2 data. 



Region p( vi CQ2- 1) y3( 12 C0 3-2) 



L1448 


2.96 


±0.42 


3.41 = 


t 0.16 


L1455 


2.86 


±0.09 


2.85 = 


t 0.30 


NGC 1333 


3.76 


±0.48 


3.52 = 


t 0.11 


B 1 


3.14 


±0.29 


3.00 = 


t0.20 


B 1 EAST 


3.16 


±0.09 


3.39 = 


t0.09 


B3 


3.36 


±0.09 


3.14 = 


t0.06 


IC 348 


2.71 


±0.42 


3.06 = 


t 0.24 
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Fig. 3. 13 CO 2-1 velocity channel maps of the Perseus region. The velocity range runs from 3 km s _1 to 1 1 km s _1 with an interval 
of 1 km s _1 which is indicated on the top of each plot. The intensities are plotted from 0.7 Kkm s _1 (~ 1 <x) to 15 Kkm s _1 . 



4.3. Velocity channel maps 

When performing the A-variance analysis not only for maps 
of integrated intensities, but for individual channel maps we 
obtain additional information on the velocity structure of the 
cloud. In the velocity channel analysis (VCA), introduced by 



Lazarian & Pogosyan (2000), the change of the spectral index 
of channel maps as a function of the channel width was used 
to simultaneously determine the scaling behavior of the density 
and the velocity fields from a single data cube of line data. Here 
we conduct such a study for the KOSMA CO data. 
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Fig. 5. A-variance analysis of the off-line channel maps. In the 
upper plot a velocity span corresponding to the integrated in- 
tensity maps is used. The two regions representing opposite ex- 
tremes in the structural behaviour, NGC 1333 and L1455, show 
about the same spectral index of the colored noise in both tran- 
sitions for small lags. In the lower plot, the influence of dif- 
ferent velocity spans, as used in the velocity channel analysis 
("Sect. l43t . is studied for L1455. The colored noise index ii no j se 
is nearly constant independent of species, transition, velocity 
range Av c h, and center velocity v cen - 



We start with the analysis of individual channel maps as 
they are provided by the channel spacing Av c h of the backends 
(cf. For all channel maps we perform the A-variance anal- 
ysis and fit power laws to the measured structure for all lags 
between the telescope beam size and the maximum scale re- 
solved by the A-variance (about 1/4 of the map size). As a re- 
sult we get the power-law index as a function of the channel 
velocity, a curve which we call index spectrum. As an example 
we show the index spectrum obtained for the 13 CO 2-1 data in 
the L 1455 region in Fig.0 The spectrum is always truncated 
at velocities where the average line temperature is lower than 
the noise rms. 

The overall structure of the index spectrum is similar to the 
line profile. The largest spectral indices are found at velocities 
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Fig. 6. A-variance spectra of integrated intensities, a) Spectra 
obtained from the CO maps and the Ay data of the region 
mapped with the KOSMA telescope, b) Spectra of integrated 
intensity maps of two 50'x 50'sub-regions: NGC 1333 and 
L 1455. Power-law fits to the data corrected for noise and beam- 
blurring are indicated as solid lines. 



close to the average line peak. This may implicate that extended 
smooth structure provides the major contribution to the overall 
emission, while the velocity tail of this structure is formed by 
small-scale features. However, the indices show an asymmetric 
behavior with respect to the blue and the red wing. The indices 
drop steeply to a noise-dominated value at the red wing, while 
the blue wing shows only a very shallow decay. Even at the 
noise limit, noticeable structure is detected in the channel maps 
there. 

For the full velocity channel analysis, the index spec- 
trum has to be computed for different velocity channel widths 

J — I 7J. 

(Lazarian & Poaosvan 2000). Thus we have binned the data to 
averages of three, five, and seven velocity channels and com- 
puted the index spectra for these binned channel maps. In Fig. [3] 
we show the results for three examples: IC 348, NGC 1333 and 
L 1455. For the sake of clarity the error bars of the index spec- 
tra were omitted in these plots. 
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Fig. 8. Comparison of the index spectra obtained for different velocity channel widths with the average line profile. The upper 
plots show the results for IC 348, the central plot NGC 1333 and the lower plot L 1455. For the left column we used the 12 CO 3-2 
data, the right column represents the 13 CO 2-1 data. The different symbols indicate the results from different velocity channel 
widths. The dashed lines represents the index of the integrated intensity maps. 



The overall structure of the index spectra is similar to Fig.0 
for all sources, transitions and channel widths. In most cases 
we find the asymmetry of a shallower blue wing relative to 
the red wing. When looking at narrow velocity channels, we 
find a dip in the centre of the index spectrum for the 12 CO 3- 
2 data of NGC 1333 and L 1455. A slight indication of such a 
dip is also present in the 12 CO 3-2 data if IC348 and in the 
13 CO 2-1 data of NGC 1333. This is due to optical depth ef- 
fects. When we check individual spectra in those regions, we 
see self-absorption in several positions. This leads to a more 



filamentary appearance of the central channel maps reflected 
by this dip in the index spectra. It is interesting to notice that 
the VCA is more sensitive to self-absorption than the average 
spectrum. 

When increasing the channel width by binning, the self- 
absorption dip is smoothed out, so that the resulting index spec- 
tra peak again close to the peak velocity of the line temperature. 
In all situations where the self absorption is negligible, the in- 
dices for the line core channels are almost independent from the 
channel width. The indices for the line integrated intensities al- 
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4 6 

velocity (kms 1 ) 

Fig. 7. Comparison of the index spectrum of the 13 CO 2-1 data 
in L 1455 with the average line profile. The index spectrum is 
created by power-law fits to the A-variance spectrum of indi- 
vidual channel maps (Av c h = 0.22 km s -1 ). The vertical error 
bars represent the uncertainty of the fit. The horizontal error 
bars indicate the velocity channel width. 



ways fall slightly below the peak indices, as they represent an 
average which is typically dominated by the line cores. 

In the red line wings, most indices remain approximately 
constant when increasing the velocity width, except for the 
largest bin width where the contribution from the core leads to 
an observable increase. In the blue wing, we find a monotonic 
growth of the spectral indices with the channel width for both 
tracers in all three regions. The additional peak at 2 km s _1 vis- 
ible in the 12 CO 3-2 data of NGC 1333 stems from a separate 
dark cloud which is also contained in the NGC 1333 map. 

Figure [9] summarizes the relation between the spectral in- 
dices and the velocity channel width. In Fig. [9^ we plot the 
average spectral index over the line as a function of the chan- 
nel width for the six data sets presented in Fig. [8] Figure |9j) 
contains the analysis when restricted to a 2 kms -1 window in 
the blue line wings. The error bars contain the standard devi- 
ation of the index variation across the line and the fit errors. 
They are necessarily large because of the systematic variation 
of in dices over the veloci t y range. In contrast to sim i lar an aly- 
sis bv lDickev et alJJ200ll) : IStanimirovic & LazarianlfeOOll) we 
find no significant systematic variation of the mean line index 
as a function of channel width (Figure |9^). In contrast to the 
average of the index spectrum we find a continuous increase 
of the spectral index with the channel width when restricting 
the analysis to the blue wing (Figure [9j>). The average index 
steepens from about 2.8 to about 3.1 in the 12 CO 3-2 maps 
and from about 2.4 to about 2.8 in the 13 CO 2-1 maps. As dis- 
cussed above we find no systematic trend in the red wings. This 
indicates that the average spectral index taken over the full line 
profile provides no measure for the velocity structure in our 
CO maps while the peculiar behaviour in the blue wings needs 
further investigation. 
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Fig. 9. Average spectral indices of the channel maps as a func- 
tion of the channel width, a) shows the average over the full line 
width, b) represents only the indices in the blue ling wings. The 
line wing components are centered at 5kms~' for NGC 1333, 
at 4^8-' for L 1455, and at 7kms^ for IC348. The error 
bars contain the fit error and the standard deviation of the in- 
dices within the considered velocity range. To avoid overlap- 
ping error bars in the plot, we have shifted the points for IC 348 
and NGC 1333 by +0.02 km s _1 relative to their actual position. 



5. Discussion 

5.1. Integrated intensity maps 

Besides the A-variance, other tools have been used to charac- 
terize interstellar cloud structure. The second-order structure 
function for an observable s(r) is S2 = (\s(r) - s(r + 6r)\ 2 ) r 
which is treat ed as a function of the absolute value of the in- 
crement \6r\ felmegreen & ScalolEooH) . IPadoan et all J2003a) 
computed the structure function of the integrated intensity map 
of 13 CO 1-0 in Perseus. A power-law fit to S2 °c Sr^ over a 
range of 0.3 to 3pc provided an index ( of 0.83. The index 
of the structure function is rel ated to the power spectral index 
by £ = p - 2 for 2 < 8 < 4 dStutzki et al.l | l998l) resulting in 
ft - 2.83. This result of Pado an et al . (2003a) agrees within the 
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error bars with the indices found by the A-variance analysis of 
the Perseus maps of integrated CO intensities over almost the 
same linear range (see Table^. 

However, the A-variance spectra of individual regions show 
significant variations of the spectral index as discussed above. 
For 13 CO 2-1, these span the range between f3 — 2.86 in L1455 
and 3.76 in NGC 1333 (Tabled). The analysis of different sub- 
sets in molecular cloud complexes thus provides additional and 
complementary information on the structure of the cloud com- 
plex. 

5.2. Velocity channel maps 

The velocity channel analysis (VCA), was used previously by 
iDickev et al.l d200ll) to study H i maps of two regions in the 4th 
Galactic quadrant. One of the regions is rich in warm H i gas, 
the ot her is rich in cool H i gas. For the warm gas, Di ckev et alJ 
find a systematic increase of the mean index with ve- 
locity channel width. The cold gas at lower latitudes behaves 
differently and shows rather constant indices of 2.7-3.1. 

The latter results resemble the outcome of the VCA of the 
12 CO and 13 CO channel maps in Perseus presented above. We 
find similar indices for the velocity integrated maps of the full 
region (cf. Tabled and the CO data show no significant vari- 
ation of the index with velocity channel width when averaged 
over all velocity bins (Figure |9^). The indices stay relatively 
constant. Since the bulk of the molecular gas traced by CO is 
even colder than cold H i gas, these results suggest a sequence 
of a reduced dependence of the spectral index of the chan- 
nel maps on the bin widths from the warm to the cold ISM. 
The constancy of the average sp ectral index could also be ex- 
plained by optical depth effects. Lazaria n & Pogos van (2004) 
have shown that absorption can lead to an effective slice broad- 
ening, which leads in extreme cases to slice indices that become 
independent from the actual channel width. 

Studying the spectral index of individual velocity bins of 
the CO data across the line profile (Fig.[8), we find that the 
power-law indices increase with the velocity channel width in 
the blue wing (Fig.|5p) while staying rather constant in the red 
win g. No corresp onding analysis was conducted for the H i data 
by D ickev et al.ld200ll) . One possibility to explain the asymme- 
try between the blue and red wings might be a shock expansion 
of the CO gas in Perseus. ISancisildl974l) found an expanding 
shell of neutral hydrogen which was created by a supernova in 
the PerOB2 association a few 10 6 years ago. At the location of 
the molecular cloud complex, the expansion is directed away 
from the Sun. Most of the associated molecular gas has been 
swept up by the shock, but pillar-like filaments have been left 
at the backside of the shock. They are visible in the channel 
maps (Fig. [3} and produce the velocity dependence seen in the 
VCA of the blue line wings, i.e. the increase of indices with 
size of the velocity bin because of a gradual increase of large 
scale contributions across the line wing. We illustrate this sce- 
nario in Fig. 1101 

The quantitative results from the velocity channel analysis 
can be interpre ted in terms of the power s pectrum of the ve- 
locity structure. Lazarian & Pogosvan I f2000 ) showed that the 



Perseus 

Molecular Cloud 



to Sun 
(350pc) 



20pc 




Fig. 10. Sketch adopted from Fig. 3 o f lSancisil ill 9741) illustrat- 
ing the spatial arrangement and motion of the Perseus cloud 
complex. The gas is swept up by a shock expansion with 
12 kms 1 . Due to the overall curvature, the line-of-sight ve- 
locity is 9 kms" 1 for IC 348, but only 7 kms" 1 for NGC1333. 
The diameter of the cloud is ~ 20 pc. Pillar-like structures are 
left at lower velocities as remainders of high-density regions 
which were not accelerated to the same velocity. 



spectrum of velocity slices as a function of the velocity channel 
width is determined by the power spectral indices of the density 
structure (3 and the velocity structure m. They obtain different 
regimes for shallow (J3 < 3) and steep (J3 > 3) density power 
spectra: 



P(k) 



P(k) oc . 




, thin slices, 
thick slices, (J3 < 3) 
very thick slices; 

thin slices, 

thick slices, (fi > 3) 
very thick slices. 



Thin slices have a velocity width less than the local veloc- 
ity dispersion at the studied scale; thick slices have a width 
larger than the velocity dispersion and very thick slices essen- 
tially correspond to the integrated maps ( Lazarian & Pogo svaiil 
2000). We can assume that a single channel of our data corre- 
sponds to thin slices as they are much narrower than any ob- 
served line width. We use j3 from the integrated maps and the 
index measured in the single channels to derive m. From the in- 
dices obtained by averaging over the full line profile (Fig.|9^), 
we obtain m values of 3.9+1.6, 3.9+1.9 and 3.8+1.7 for 12 CO 
3-2 in IC348, NGC 1333 and L 1455, respectively; while m 
is 3.9+2.0, 3.5+2.5 and 4.1 + 1.8 for 13 CO 2-1 in the same 
three regions. In the blue wings (Fig. |9j>), we obtain m val- 
ues of 3.2+1.6, 3.4+1.4 and 3.6+1.5 for 12 CO 3-2; while m is 
4.3± 1 .9, 4.3+ 1 .7 and 4. 1 ± 1 .5 for I3 CO 2-1 in the three regions. 
These values have large error bars, so that they are not directly 
suited to discriminate between different turbulence models. At 
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least, we find that all values are consisten t with Kolmogorov 
turbulence that gives m ~ 3.7 ( Kol mo gorovl 194-iT) . 

6. Summary 

1. We present KOSMA maps of the I3 CO 2-1 and I2 CO 3-2 
emission of the Perseus molecular cloud covering 7.1 deg 2 . 
These data are combined with FCRAO maps of integrated 
12 CO and 13 CO 1-0 intensities and with a 2MASS map of 
optical extinctions. 

2. To characterize the cloud density structure, we applied the 

A-variance analysis to integrated intensity maps. The A- 

variance spectra of the overall region follow a power law 

with an index of B — 2.9 - 3.0 for scales between 0.2 

r I 1 

and 3 pc. This agrees with results obtained by Padoan et al. 

(2003a) studying structure functions of a 13 CO 1-0 map of 

Perseus. 

3. We also applied the A-variance method to seven sub- 
regions of Perseus. The resulting power spectral indices 
vary significantly between the individual regions. The ac- 
tive star-forming region NGC 1333 shows high spectral in- 
dices (J3 = 3.5 - 3.8) while the dark cloud L 1455 shows 
low indices of 2.9 in both transitions. 

4. Additional information is obtained from the A-variance 
spectra of individual velocity channel maps. They are very 
sensitive to optical depth effects, indicating self-absorption 
in the densest regions. The asymmetry of the channel map 
indices relative to the line centrum is a hint towards a pecu- 
liar velocity structure of the Perseus cloud complex. 

5. When analyzing the spectral indices as a function of the ve- 
locity channel width we find almost constant indices when 
averaging over the total line profile. A continuous increase 
of the index with varying velocity channel width is, how- 
ever, observed in the blue wings. This behavior can be ex- 
plained by a shock running through the region creating a 
filamentary structure preferentially at low velocities. 

We find that the comparison of the structural properties 
for entire surveys and sub-sets, as well as the velocity channel 
analysis (VCA), provide additional, significant characteristics 
of the ISM in observed CO spectral line maps. These quantities 
are useful for a comparison of the structure observed in differ- 
ent clouds, possibly providing a diagnostic tool to characterize 
the star-formation activity and providing additional constraints 
for numerical simulations of the ISM structure. 
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